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ABSTRACT 



A method and apparatus are provided that allows initial 
beamforming weights to be determined before timing or 
delay and frequency offset are known. According to one 
aspect of the present invention, the invention includes 
receiving a burst from a terminal at a set of spatial diversity 
antennas, the burst including a known training sequence and 
storing measurements of the received burst from each spatial 
diversity antenna corresponding to receive times within a 
selected analysis window. The invention further includes 
determining weights for each diversity antenna for at least a 
portion of the training sequence of the burst before timing 
for the burst has been determined, applying the determined 
weights to the stored measurements and determining timing 
for the burst based on the stored measurements to which the 
determined weights have been applied. 

31 Claims, 7 Drawing Sheets 



V V V v 



[■"MQ3 



DUPLEXER 



4iRF 
RECEIVER 
MODULES 



-205 



ADC 

IE 



y-209 



DOWN 
CONVERTER 



(4x4 



-107 



RF TRANSMIT 
MODULES 



RF/TWING 
CONTROLLER 



^233 



TIMESLOT 
PROCESSORS 



-245 



TRANSMIT 
CONTROLLER/ 
W00ULAT0R 



-237 



HOST 
OSP 



^231 



HIGHER LEVEL PROCESSING 



06/14/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 18,2003 Sheet 1 of 7 US 6,650,881 Bl 



FIG. 1 



V V W V 



-103 




-205 



4xRF 
RECEIVER 
MODULES 



I 



ADC 

DOWN 
CONVERTER 

{4x4 

TIMESLOT 
PROCESSORS 



-209 



"_^-213 



-217 



DUPLEXER 

RX TX 



RF/TIMING 
CONTROLLER 



HOST 
DSP 



-107 



-245 



RF TRANSMIT 
MODULES 



-233 



TRANSMIT 
CONTROLLER/ 
MODULATOR 



y-237 



231 



HIGHER LEVEL PROCESSING 



06/14/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 18, 2003 Sheet 2 of 7 US 6,650,881 Bl 



FIG. 2 



V 



61 

A 



.45 



DUPLEXER 



if 



46 



I 



47 



TRANSMITTER 



59- 



60 



62 



68- 



57- 



MODULATOR 



63 

A 



CPU 



58 



•65 



64- 



MICROPHONE 



56- 




50 

A 



RECEIVER 



—48 



s 

52 



49 



DEMODULATOR 



-51 



67 



SPEAKER 



■66 



■54 




— 53 



06/14/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 18, 2003 Sheet 3 of 7 US 6,650,881 Bl 



10[is 



272 ms 



256 us 



16|is 



S12ms 



10ms 14ms 




CORRECTION CORRECTION S 
TRAINING Wfi$ TRAINING $ 



«•••*•»« 



1090ms 



■: BROADCAST 

(BCH) 
:■ INFORMATION 



I I 



I 



FIG. 3 



10^ 



240 



164ns 10 |is 106jis 15ps 




fig. 4 



06/14/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 18, 2003 Sheet 4 of 7 US 6,650,881 Bl 



520 



r 

Roundtrip Delay Margin Window 



r 



530 



r 



505 



Beamfbrmer Analysis | Marker t 
Window i I 























Pi 


P2 


P3 


Pa 


Ps 


Pe 


P7 


Pa 


P9 


:< : :M:-x 





















[X t , Xj, X2, X 3 , X4, X 5 , Xg, X 7 , X B , X 9 . X 10 . X 11t X12J 



510 



J 



[-Xl -X,, -X 2 , -X 3 , -X4, -X 5 , -X 6 , -X 7 , -Xe. -Xg, -x 10 , -xn, -x 12 ) 



FIG. 5 



06/14/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 18, 2003 Sheet 5 of 7 



US 6,650,881 Bl 



FIG. 6 



RECEIVE THE CR BURST AT THE ANTENNA ARRAY 



STORE THE MEASUREMENTS FROM THE 
MULTIPLE ANTENNAS 



FIND A COARSE RELATIVE TIMING HYPOTHESIS 



PERFORM A FOLLOW-UP RELATIVE TIMING SEARCH 



I 



DETERMINE A WEIGHT FACTOR 



CONVERT THE MEASUREMENTS FROM EACH ANTENNA 
CHANNEL INTO A SINGLE CHANNEL 



I 



DETERMINE THE FINE TIMING 



DETERMINE THE FREQUENCY OFFSET 



± 



DETERMINE A MORE ACCURATE WEIGHT VECTOR USING 
DATA IN THE BEAMFORMING ANALYSIS WINDOW 



CONVERT THE MEASUREMENTS FROM EACH ANTENNA 
CHANNEL INTO A SINGLE CHANNEL 



DETERMINE ABSOLUTE TIMING 



DETERMINE A MORE ACCURATE WEIGHT VECTOR USING 
THE FULL TRAINING SEQUENCE 



CONVERT THE MEASUREMENTS FROM EACH ANTENNA 
CHANNEL INTO A SINGLE CHANNEL 



DEMODULATE THE CR BURST USING THE SINGLE CHANNEL 



-602 
-604 

-606 
-608 
-610 
-612 

-614 
-616 
-618 

f—620 

-622 
-624 

-626 

-628 



06/14/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 18,2003 Sheet 6 of 7 



US 6,650,881 Bl 



FIG. 7 



IS ERROR 

t 

I 

^-702 



X 



0 2/3 



coreSnapPoint 



706 



4/3 



C 



704 



(— ^ TIME DELAY, X 



11 1/3 



BEST HYPOTHESIS 
FOR COARSE SEARCH 



06/14/2004, EAST Version: 1.4.1 



U.S. Patent Nov. 18, 2003 Sheet 7 of 7 US 6,650,881 Bl 



FIG. 8 



IS ERROR f(t x ) 



f(11 1/3) 



t 



802 



806 

ill! 

I I 



I 



804 



-2/3 



I 

0 

I 

coreSnapPoint 



2/3 i 1 



! 4/3J 



X BEST HYPOTHESIS 

FROM COARSE SEARCH 



BEST HYPOTHESIS 
FROM FOLLOW-UP 



06/14/2004, EAST Version: 1.4.1 



US 6,6f 

1 

CALCULATING SPATIAL WEIGHTS IN A 
RADIO COMMUNICATIONS SYSTEM 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention applies to the field of spatial 
diversity radio communications systems and, in particular, 
to determining spatial weights for a received burst for which 
timing are unknown. 

2. Description of the Prior Art 

Mobile radio communications systems such as cellular 
voice and data radio systems typically have several base 
stations in different locations available for use by mobile or 
fixed user terminals, such as cellular telephones or wireless 
web devices. Each base station typically is assigned a set of 
frequencies or channels to use for communications with the 
user terminals. The channels are different from those of 
neighboring base stations in order to avoid interference 
between neighboring base stations. As a result, the user 
terminals can easily distinguish the transmissions received 
from one base station from the signals received from 
another. In addition, each base station can act independently 
in allocating and using the channel resources assigned to it. 

Such radio communications systems typically include a 
broadcast channel (BCH). The BCH is broadcast to aU user 
terminals whether they are registered on the network or not 
and informs the user terminals about the network. In order 
to access the network, a user terminal normally tunes to and 
listens to the BCH before accessing the network. It will then 
use the information in the BCH to request access to the 
network. Such a request typically results in an exchange of 
information about the network using separate control and 
access channels and ends in the user terminal receiving an 
assignment to a particular base station. 

While frequency and timing oflset can sometimes be 
determined by a user terminal based on the BCH, the initial 
request for access is typically transmitted by the user ter- 
minal without any knowledge of its distance to the base 
station, its timing delay or the direction to the base station. 
In a spatial diversity multiple access system, the base station 
can enhance the capacity of the system by determining the 
position and range to the user terminal as well as any other 
spatial parameters. The timing uncertainty of the arrival time 
of such request messages are proportional to the round trip 
delay encountered by messages traveling between the base 
station and the mobile terminal. For systems with a high 
coverage area per base station, this range and therefore the 
delay uncertainty may be very large. For example, a range 
of fifteen km results in a roundtrip delay time of around 100 
microseconds. 

In order to accurately resolve the access request and 
determine spatial parameters the timing or delay, frequency 
offset and spatial signature of the user terminal's message 
are all desired. Typically, it is preferred to accurately deter- 
mine the timing, frequency oflset and delay before the 
spatial signature is determined. This can create great 
demands on the processing resources of a base station and 
increase the amount of time required to generate a reply to 
the access request. 

BRIEF SUMMARY OF THE INVENTION 

A method and apparatus are provided that allows initial 
beamforming weights to be determined before timing or 
delay and frequency offset are known. According to one 
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aspect of the present invention, the invention includes 
receiving a burst from a terminal at a set of spatial diversity 
antennas, the burst including a known training sequence and 
storing measurements of the received burst from each spatial 

5 diversity antenna corresponding to receive times within a 
selected analysis window. The invention further includes 
determining weights for each diversity antenna for at least a 
portion of the training sequence of the burst before timing 
for the burst has been determined, applying the determined 

10 weights to the stored measurements and determining timing 
for the burst based on the stored measurements to which the 
determined weights have been applied. 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 The present invention is illustrated by way of example, 
and not by way of limitation, in the figures of the accom- 
panying drawings in which like reference numerals refer to 
similar elements and in which: 
2Q FIG. 1 is a simplified block diagram of a base station on 
which an embodiment of the invention can be implemented; 

FIG. 2 is a block diagram of a remote terminal on which 
an embodiment of the invention can be implemented; 
FIG. 3 is a diagram of an exemplary broadcast channel 
25 BCH burst; 

FIG. 4 is a diagram of an exemplary configuration request 
CR burst; 

FIG. 5 is a diagram of an exemplary training sequence in 
accordance with the present invention; 
30 FIG. 6 is a flow chart of an implementation of aspects of 
the present invention; 

FIG. 7 is a graph of a search for a least squares error 
relative timing; and 
35 FIG. 8 is a graph of a follow-up search for a least squares 
error relative timing. 

DETAILED DESCRIPTION OF THE 
INVENTION 

40 Base Station Structure 

The present invention relates to wireless communication 
systems and may be a fixed-access or mobile-access wireless 
network using spatial division multiple access (SDMA) 

45 technology in combination with multiple access systems, 
such as time division multiple access (TDMA), frequency 
division multiple access (FDMA) and code division multiple 
access (CDMA). Multiple access can be combined with 
frequency division duplexing (FDD) or time division 

50 duplexing (TDD). FIG. 1 shows an example of a base station 
of a wireless communications system or network suitable for 
implementing the present invention. The system or network 
includes a number of subscriber stations, also referred to as 
remote terminals or user terminals, such as that shown in 

55 FIG. 2. The base station may be connected to a wide area 
network (WAN) through its host DSP 231 for providing any 
required data services and connections external to the imme- 
diate wireless system. To support spatial diversity, a plurality 
of antennas 103 is used, for example four antennas, although 

60 other numbers of antennas may be selected. 

A set of spatial multiplexing weights for each subscriber 
station are applied to the respective modulated signals to 
produce spatially multiplexed signals to be transmitted by 
the bank of four antennas. The host DSP 231 produces and 

65 maintains spatial signatures for each subscriber station for 
each conventional channel and calculates spatial multiplex- 
ing and demultiplexing weights using received signal mea- 
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surements. In this manner, the signals from the current active 
subscriber stations, some of which may be active on the 
same conventional channel, are separated and interference 
and noise suppressed. When communicating from the base 
station to the subscriber stations, an optimized multi-lobe 
antenna radiation pattern tailored to the current active sub- 
scriber station connections and interference situation is 
created. Suitable smart antenna technologies for achieving 
such a spatially directed beam are described, for example, in 
U.S. Pat. Nos. 5,828,658, issued Oct. 27, 1998 to Ottersten 
et al. and 5,642,353, issued Jun. 24, 1997 to Roy, III et al. 

The outputs of the antennas are connected to a duplexer 
switch 107, which in this TDD system is a time switch. Two 
possible implementations of switch 107 are as a frequency 
duplexer in a frequency division duplex (FDD) system, and 
as a time switch in a time division duplex (TDD) system. 
When receiving, the antenna outputs are connected via 
switch 107 to a receiver 205, and are mixed down in analog 
by RF receiver ("RX") modules 205 from the carrier fre- 
quency to an FM intermediate frequency ("IF"). This signal 
then is digitized (sampled) by analog to digital converters 
("ADCs") 209. Only the real part of the signal is sampled. 
Final down -converting to baseband is carried out digitally. 
Digital filters can be used to implement the down-converting 
and the digital filtering, the latter using finite impulse 
response (FIR) filtering techniques. This is shown as block 
213. The invention can be adapted to suit a wide variety of 
RF and IF carrier frequencies and bands. 

There are, in the present example, four down-converted 
outputs from each antenna's digital filter device 213, one per 
receive timcslot. The particular number of timeslots can be 
varied to suit network needs. While the present example uses 
four uplink and four downlink timeslots for each TDD 
frame, desirable results have also been achieved with three 
timeslots for the uplink and downlink in each frame. For 
each of the four receive timeslots, the four down-converted 
outputs from the four antennas are fed to a digital signal 
processor (DSP) device 217 (hereinafter "timeslot 
processor") for further processing, including calibration, 
according to one aspect of this invention. Four Motorola 
DSP56303 DSPs can be used as timeslot processors, one per 
receive timeslot. The timeslot processors 217 monitor the 
received signal power and estimate the frequency offset and 
time alignment. They also determine smart antenna weights 
for each antenna element. These are used in the spatial 
diversity multiple access scheme to determine a signal from 
a particular remote user and to demodulate the determined 
signal. 

The output of the timeslot processors 217 is demodulated 
burst data for each of the four receive timeslots. This data is 
sent to the host DSP processor 231 whose main function is 
to control all elements of the system and interface with the 
higher level processing, which is the processing which deals 
with what signals are required for communications in all the 
different control and service communication channels 
defined in the system's communication protocol. The host 
DSP 231 can be a Motorola DSP56303. In addition, timeslot 
processors send the determined receive weights for each 
user terminal to the host DSP 231. The host DSP 231 
maintains state and timing information, receives uplink burst 
data from the timeslot processors 217, and programs the 
timeslot processors 217. In addition it decrypts, 
descrambles, checks error correcting code, and deconstructs 
bursts of the uplink signals, then formats the uplink signals 
to be sent for higher level processing in other parts of the 
base station. With respect to the other parts of the base 
station it formats service data and traffic data for further 
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higher processing in the base station, receives downlink 
messages and traffic data from the other parts of the base 
station, processes the downlink bursts and formats and sends 
the downlink bursts to a transmit controllermodulator, 
S shown as 237. The host DSP also manages programming of 
other components of the base station including the transmit 
controllermodulator 237 and the RF timing controller shown 
as 233. 

The RF timing controller 233 interfaces with the RF 

1° system, shown as block 245 and also produces a number of 
timing signals that are used by both the RF system and the 
modem. The RF controller 233 reads and transmits power 
monitoring and control values, controls the duplexer 107 
and receives timing parameters and other settings for each 

15 burst from the host DSP 231. 

The transmit controllermodulator 237, receives transmit 
data from the host DSP 231, four symbols at a time. The 
transmit controller uses this data to produce analog IF 
outputs which are sent to the RF transmitter (TX) modules 

20 245. Specifically, the received data bits are converted into a 
complex modulated signal, up-converted to an IF frequency, 
4-times over-sampled, multiplied by transmit weights 
obtained from host DSP 231, and converted via digital to 
analog converters ("DACs") which are part of transmit 

25 controllermodulator 237 to analog transmit wavefonnns. 
The analog waveforms are sent to the transmit modules 245. 

The transmit modules 245 up-convert the signals to the 
transmission frequency and amplify the signals. The ampli- 

30 fied transmission signal outputs are sent to antennas 103 via 
the duplexertime switch 107. 

User Terminal Structure 

FIG. 2 depicts an example component arrangement in a 

35 remote terminal that provides data or voice communication. 
The remote terminal's antenna 45 is connected to a duplexer 
46 to permit antenna 45 to be used for both transmission and 
reception. The antenna can be omni-directional or direc- 
tional. For optimal performance, the antenna can be made up 

40 of multiple elements and employ spatial processing as 
discussed above for the base station. In an alternate 
embodiment, separate receive and transmit antennas are 
used eliminating the need for the duplexer 46. In another 
alternate embodiment, where time division diversity is used, 

4 5 a transmitreceive (TR) switch can be used instead of a 
duplexer as is well-known in the art. The duplexer output 47 
serves as input to a receiver 48. The receiver 48 produces a 
down-converted signal 49 which is the input to a demodu- 
lator 51. A demodulated received sound or voice signal 67 

50 is input to a speaker 66. 

The remote terminal has a corresponding transmit chain in 
which data or voice to be transmitted is modulated in a 
modulator 57. The modulated signal to be transmitted 59, 
output by the modulator 57, is up-converted and amplified 

55 by a transmitter 60, producing a transmitter output signal 61. 
The transmitter output 61 is then input to the duplexer 46 for 
transmission by the antenna 45. 

The demodulated received data 52 is supplied to a remote 
terminal central processing unit 68 (CPU) as is received data 

60 before demodulation 50. The remote terminal CPU 68 can 
be implemented with a standard DSP (digital signal 
processor) device such as a Motorola series 56300 DSP. This 
DSP can also perform the functions of the demodulator 51 
and the modulator 57. The remote terminal CPU 68 controls 

65 the receiver through fine 63, the transmitter through line 62, 
the demodulator through line 52 and the modulator through 
line 58. It also communicates with a keyboard 53 through 
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line 54 and a display 56 through line 55. A microphone 64 Each broadcast message is mapped into a broadcast burst 

and speaker 66 are connected through the modulator 57 and with the information shown below in Table 2. 
the demodulator 51 through lines 65 and 66, respectively for 

a voice communications remote terminal. In another TABLE 2 

embodiment, the microphone and speaker are also in direct 5 
communication with the CPU to provide voice or data 
communications. 

The remote terminal's voice signal to be transmitted 65 
from the microphone 64 is input to a modulator 57. Traffic 

and control data to be transmitted 58 is supplied by the 1Q 
remote terminal's CPU 68. Control data 58 is transmitted to 
base stations during registration, session initiation and ter- 
mination as well as during the session as described in greater 
detail below. 

In an alternate embodiment, the speaker 66, and the ^ BStxPwr is the effective isotropic radiated power of the 
microphone 64 are replaced or augmented by digital inteT- broadcast message. This number indicates the power trans- 
faces well-known in the art that allow data to be transmitted mitted by the base station taking into account the number of 
to and from an external data processing device (for example, amplifiers and diversity antennas available at the base sta- 
a computer). Io one embodiment, the remote terminal's CPU ti 0 n. 

is coupled to a standard digital interface such as a PCMCIA BSCC ^ the Dase stau0Q coi or co de, used by the user 

interface to an external computer and the display keyboard, M terminal t0 sc i ect traimng data for uplink bursts and to 

microphone and speaker are a part of the external computer. broadcasts of different base stations. 

The remote terminal s CPU 68 communicates with these __f , . . - . 4 f , 

components through the digital interface and the external fc B k Sload 15 ™ ^ dlcatlon of te amount of ™ d ca P acit * 

computer's controller. For data only communications, the thc basc statl0n has - 

microphone and speaker can be deleted. For voice only 25 ln one embodiment, the network is designed to take 
communications, the keyboard and display can be deleted. maximal advantage of spatial division multiple access tech- 
, m «jv nologies and particularly smart antenna array signal pro- 
Broadcast Channel (BCH) ce&sing To heJp maintain reliable spatial chan nels in an 

In one embodiment, the system of the present invention is extremely dense frequency reuse pattern, the network uses 

initiated for each user terminal or remote terminal from the 3Q t | me division duplex TDM A where uplink and downlink 

broadcast channel BCH which is transmitted as a burst from transmissions are always on the same frequency. In addition, 

the base station to all potential user terminals. The BCH De cause many user terminals are single antenna and transmit 

burst, unlike the traffic channel bursts, is transmitted in all ^ fcccivc omnidirC ctionaUy, except for the BCH, an 

directions where user terminals may be, typically omnidi- {h± ^ [& ^ received 5efore a doW nlink burst 

rectionally but the specific beam pattern will depend on the ^ tQ bc ^ ^ aUows downlink bursts to bc more 

network. Anexample of a broadcast burst structure is shown accuratel atiaU direcled ^ link training sequence is 

in FIG. 3. The BCH communicates enough basic mforma- cmbcdde ^ i^ery' uplink burst to allow for moderately fast 

tion to enable a subsequent exchange of configuration f h ^ d ite dec orrelation of the spatial 

request CR and configuration message CM between the base ^ ith ^ ency . 

station and the user terminal. The BCH also provides general cc be Qnc 

frequency offset and timing update information to all user ^ nni \ q Z n ™ w Vknown in the art. In one 

terminas. embodiment, the parameters of the frequency hopping 

Table 1, below summaries the content of an example of scheme afe initiaUy unknQwn tQ ^ uger TO& 

a BCH burst, as shown in FIG. 3. maximizes the flexibility of the network and increases the 

45 flexibility of the user terminal. As explained below, the 

TABLE 1 frequency hopping parameters are transmitted to the user in 

Duration Contents " the CM burst. 

— ^ j n one embodiment, the BCH channel is shared by all base 

lO/aec ramp -up stations in the wireless communication system. Using the 7 

%Z » bit BSCC, up to 128 base stations can be accommodated. 

16 fiscc broadcast preamble r„ r* . . . , r e The BCH is part of a time division duplex channel with a 

512 //sec information symbols h\, h'* . . . , b'a* repeating frame. The channel that includes BCH is a single 

10 //sec ramp -down RF carrier frequency used for uplink and dowDlink. For high 

14 /see mter-burst guard tune ^ environments Qf (qt iDcreased robustness, the BCH 

55 can hop frequencies according to a predetermined scheme or 

The frequency and timing correction training symbols can be repeated on several different frequencies. The repeating 

be set according to any one of many appro aches well-known frame includes the downlink BCH for each base station, 

in the art. They can also be combined, exchanged with a labeled BS1 etc. as shown in Table 3 below. The next frame 

synchronization sequence or eliminated. includes the uplink Configuration Request CR, labeled CR1 

The broadcast information symbols are constructed from 60 etc. and downlink Configuration Message CM, labeled CM1 

a 15-bit broadcast message which is modulated and coded etc. Each frame also includes a number of reserved slots, 

into a 256 bit sequence. The number of symbols as well as shown as empty boxes below. These slots can be used for 

the structure and sequence of transmitted bits can be varied data traffic, if the broadcast channel is also used for traffic, 

to suit a wide variety of applications. The broadcast channel for other control messages or reserved to reduce interference 

information symbols provide the information needed for a 65 on other channels in the network. In one embodiment, the 

user terminal to request a configuration message from the other traffic channels hop frequencies around and through 

base station. the BCH. The frames are repeated for each respective base 
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station 1 to 128 to build a superframe as discussed in more 
detail below. After the last CM, CM128, the superframe 
repeats and begins again with the next superframe and the 
BCH for base station 1. 



TABLE 3 





Uplink 


Downlink 


Superframe 
1 


Frame 1 








BSl 






Frame 2 


CR1 






CM1 






Frame 3 








BS2 






Frame 4 


CR2 






CM2 












Frame 255 








BS128 






Frame 256 


CR128 






CM 128 






Superframe 
2 


Frame 1 








BSl 






Frame 2 


CfU 






CM1 













In another embodiment, the BCH is on its own channel 
and CR and CM are on a separate control channel. 
Alternately, one BCH can be provided on a constant fre- 
quency and a secondary BCH can be provided on another 
channel with hopping frequency. The hopping channel is 
described in the CM. 

Registration 

In one embodiment of the present invention, a session 
begins when a user terminal seeks to register with a base 
station. The user terminal does this without knowing its 
relative position to the best station in terms of direction and 
distance (or range) to the base station (BS). Therefore, when 
the user terminal (UT or remote terminal or subscriber 
station) requests registration using a configuration request 
(CR) burst, the base station receives the transmission from 
the user terminal with a large timing uncertainty and 
unknown spatial parameters or weights when compared to 
traffic bursts that already utilize timing advance directives 
from the base station. 

Before registration, several user terminals may transmit 
CR bursts during the same time slot to register with the same 
base station. CR bursts addressed to different base stations 
may also be received. The base station resolves these 
requests using spatial processing and multiple antennas. By 
combining the multiple antenna measurements 
(beamforming) the base station minimizes the interference 
between signals and maximizes the signal-to-noise ratio 
(SNR) of each received burst. 

Beamforming can be performed in a variety of ways. In 
one embodiment, a training sequence is used with a least- 
squares cost function to determine the weights for a beam- 
former. This allows signals to be classified as desired or 
undesired. The training sequence -based approach uses an 
estimation of the timing and frequency offset of the CR 
burst. This estimation task can be performed by determining 
the least-squares cost for each timing and frequency hypoth- 
esis. This demands great computational resources. Using a 
special sequence design and a beamforming algorithm that 
does not require a search over the entire timing uncertainty 
range alleviates the computational load of the full search, In 
one embodiment, since CR bursts have a high timing 
inaccuracy, a periodic training sequence is dedicated for CR 
bursts alone. 



Configuration Request Burst Structure 

A configuration request (CR) burst is transmitted by a user 
terminal (user terminal) in order to initiate communications 
or registration with a base station (base station). It is 

S transmitted after gathering information about the system by 
listening to (multiple) broadcast channel (BCH) bursts. The 
CR burst, is the first communication from a user terminal to 
the base station and therefore, the user terminal does not 
have any information about its range to its selected base 

10 station. Accordingly timing, range, and spatial processing 
weights among other things are all unknown to the base 
station. An example of a CR burst is shown in FIG. 4. 

The configuration request burst is composed of several 
fields as shown in FIG. 4 which are listed in Table 4. The 

35 durations are described in terms of microseconds. In one 
embodiment, a symbol period is 2 microseconds. 



25 



TABLE 4 


Configuration Request (CR) Burst fields 


Duration 


Contents 


10 face 


ramp-up 


240 //sec 


training symbols 


1 64 //sec 


information symbols 


10 //sec 


ramp - down 


106 //sec 


extra guard time 


15 //sec 


inter-burst guard time 



The training symbols are allocated 240 microseconds in 
30 order to allow the signal to be accurately received and 
demodulated before a user terminal has registered and has 
received any knowledge of the system. The training symbols 
are discussed in more detail below. 
The 82 information symbols are constructed from the 
35 configuration request message using, for example, forward 
error correction coding. In the present embodiment, the CR 
burst is modulated using ji/2-BPSK modulation in order to 
decrease the peak- to -average ratio of the transmitted wave- 
form. 

The information symbols of the present CR burst are 
40 mapped out as shown in Table 5, below. Any of the items 
listed below can be deleted and transmitted later during the 
registration cycle or not at all based on the needs of the 
system. CR is scrambled by a function of BSCC ensuring 
that even if there is some interference from CRs sent to 
45 nearby base stations, the demodulation capture effect of the 
BSCC works out any collisions. In one embodiment, the 
scrambling is performed by taking the encoded bit sequence 
and exclusive OR'ing it with the output of a linear feedback 



50 



shift register. 






TABLE 5 


Configuration Request Message 


Field 


#of Bits 


identity 


8 


utClass 


4 


txPwr 


5 


Total 


17 



identity is a set of unique random bits for each user 
terminal that differentiate simultaneous messages from 
multiple user terminals. Because of the randomness and 
65 large number of bits, it is unlikely that two user 
terminals will select the same identity code at the same 
time. 
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utClass identifies user terminal capabilities (highest the training sequence. The period of the core sequence is a 

modulation class, frequency hopping capabilities, etc.) compromise between the duration of the search range and 

This sequence identifies the type of user terminal that the possibility of finding sequences of the specified period 

sent the CR. A palmtop digital assistant might have with good autocorrelation and cross correlation properties, 

different capabilities than a desktop computer with a 5 In the present embodiment, constraints on the periodicity, 

fixed dedicated antenna. With utClass, the different mean, autocorrelation and cross correlation are used to 

capabilities can be distinguished. define the core sequences of the training sequence for CR 

txPwr represents the power used by the user terminal to hurst acquisition. 

transmit the Configuration Request burst. For example, Periodicity: x(k+P)=x(k), in which x(k) is the periodic 

user terminal power=(2txPwr-30) dBm. 10 p ar t of a n/2-BPSK training sequence and P is the period of 

In the present embodiment, the CR burst includes 106 the sequence. In the present example, P=12 symbols. 

microseconds of extra guard time. As mentioned above, the Mean: the mean of the periodic sequence is constrained to 

CR burst is transmitted by a user terrninal without any be zero: 

knowledge of the time required for the burst to travel to the 

base station. In the present embodiment, the desired maxi- 15 ^ p 

mum range from a user terminal to a base station is 15 km. x(k) = 0. 

The delay of the burst's reception at the base station includes *=' 

the delay in receiving the BCH from the base station. A user 

terminal at range R km receives the BCH burst(s) R/c . . c , f , . , x 

. ~ ., & . ... , , , on Autocorrelation: The out-of phase (symbol spaced) auto - 

seconds after they are transmitted by the base station (c=3x 20 r w r / 

10° m/sec). When the user terminal decides to transmit a CR ™ TT ^°? » bounded b * 1/3 > where T ,S the P hase oflket of 
burst, its signal wiU be delayed by another R/c seconds, so ™t-of-phase sequence: 

that the delay of the CR burst at the base station relative to 
base station time will be 2R/c seconds. To accommodate 
delays created by a maximum range of 15 km, the guard 25 

period at the end of the CR burst time slot should be at least p " * 1 /3 * T * 0 

100 microseconds. In addition, in the present embodiment, 2, W*)l 

the user terminal randomly delays its CR burst transmission 
from 0 to 9 symbols in order to reduce collisions with other 

user terminals that are sending CR bursts in the same slot 30 Cross correlation: Symbol spaced cross correlation 
and that may be at about the same range from the base between two candidate periodic sequences x(k) and y(k) is 
station. Therefore, a CR burst from a distant user terminal bounded by 1/3: 
may arrive at the base station as much as 118 microseconds 
late as compared to a user terminal that is very close to the 
base station and that does not delay its CR burst transmis- 35 
sion. This accounts for the 106 microseconds of extra guard 



time that is allotted to the CR burst and the extra margin / £ ^ k ^z / £ 

supplied by the 15 microseconds of inter-burst guard time. V *=i V *=i 
With the inter-burst guard time 3 microseconds of extra 

margin are allotted before the beginning of the next slot. 4 ° DDCV * t iU 

* < j • . 1 • 1 j aau- 1 1c At least two ti/2-BPSK sequences conform to these con- 

The present embodiment also includes an additional 15 , . t u ui #u ♦ 

F , f . , . . , ,. u . , , . strain ts. More sequences may be possible if the constraints 

microseconds of interburst guard time. J nis helps to confine , - n ^ ^ r 

the RF energy from the CR burst to the time-slot allocated arc f™*' Thcsc s ^™ a ™ ™ us f cd as * c corc H 12 

for the CR burst in order not to create adjacent slot inter- ^ mbo{ sc « uence *J 15 re P eat f to ^ orm the ^ no6lC 

ference. The interburst guard time provides a total guard « 'raining sequence as discussed below. They are: 

time of 121 microseconds to reduce the likelihood that a CR s^flj,! j,l,-j,-ld,-l,-j,-l,-jj 

burst will extend into the next time slot. s 2 o[l,j,l,-j,-l,-j,l 

The CR burst is sent on the control carrier, as an example, The bounds on auto and cross correlations help to make 

exactly 2265 /(see after receipt of a downlink BCH burst. In delayed versions of these sequences to appear partially 

this way, an otherwise uninitialized user terminal can send 5 0 uncorrelated to an LS beamformer which resolves them. 

CR without any knowledge of the frequency hopping The two periodic sequences can be shared among base 

sequence parameters. As discussed above, the CR burst is stations. In one embodiment, depending on its color code, a 

shorter than a standard uplink time-slot to allow for base station will accept one of the two core training 

unknown time-of-flight from the user terminal to the base sequences but not both. The user terminal will be able to 

station and typically arrives late in the uplink time-slot 55 transmit the correct accepted sequence after learning the 

receive window. base station color code from the broadcast (BCH) burst. This 

enables a reuse pattern of 2. For reuse patterns greater than 

Training Sequences 2 , some of the constraints on the sequences can be relaxed 

In the present embodiment, known periodic training to yield a larger set of possible periodic sequences to 

sequences are used for the training of base station beam- 60 distribute among the base stations. 

formers in order to extract the desired CR bursts in the In an alternate embodiment, the user terminal selects 

presence of interference. The training sequence includes which of the core sequences it will transmit. This can aid a 

several periodic repetitions of a core sequence followed by base station in resolving two interfering CR bursts. The core 

a marker sequence. Alternatively, the marker sequence may sequence can be selected based on a serial number, product 

precede the repetition of the core sequence. By repeating the 65 number, ID number or other stored number of the user 

core sequence, the search range for a least squares (LS) terminal. For, example if the two sequences s 1 and Sj above 

beamformer can be reduced to a single repetition period of are available, the user terminal can read its serial number 
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register. If the serial number is odd it selects sequence si and of x; P1,P2, . . . , P9 are followed by a marker sequence M 

if the serial number is even it selects S 2 . In another 505. FIG. 5 also can be viewed from the perspective of base 

embodiment, the user terminal can generate a random num- station timing. In one embodiment, FIG. 5 shows windows 

ber from 1 to 2. If there are more than two sequences a that the base station uses to resolve received training 

greater range of random numbers can be generated. The s sequences in a CR burst. In this context, the sequences PI to 

random number is then mapped into a list of sequences in M would line up with the windows as shown only if the 

order to select the periodic sequence that is associated with sequences were received from a user terminal that was in the 

the generated random number. same location as the base station, or in other words, if there 

were no propagation delay between the base station and the 
Transmitting CR Burst and Random Delay 10 user terminal. In this paradigm, a round trip delay margin 
In one embodiment, the base station transmits broadcast window 520 contains the sequences PI through P5 for a total 
channel (BCH) bursts continuously over a pre-established of 60 symbols or 120 microseconds. This accommodates a 
known channel. To begin a session with a desired base combined delay of 50 symbols due to the round trip propa- 
gation, an unregistered and unrecognized user terminal in an g ation dela V between the user lerrainal and the base j s ^ on 
unknown location listens to the BCH bursts. From the BCH, 15 and 10 symbols of user terminal selected random delay. A 
it determines the base station with which wants to establish beamformer analysis window 530 contains core sequences 
communications and selects an appropriate training P6 through P9 and so is 48 symbols wide. This window will 
sequence for that base station. It also selects a random capture 4 periods of x even if the CR burst is delayed by up 
number D between 0 and 9, and transmits its CR burst after t0 60 symbols. 

delaying it for D symbols. The number for the delay can be 20 Using the observations that fall into the beamformer 

selected as described above with respect to the periodic analysis window, a base station can estimate the delay 

training sequence. The user terminal can select an internal encountered by the CR burst modulo-12 symbols since each 

register such as a Serial No. or Product ID No., that contains period PI through P9 is identical and 12 symbols long. In 

a number greater than 9 and use this as a pointer to select the other words, the measurement in the beamformer analysis 

number D. Alternatively, the user terminal can use a random 25 window will be the same for the delays % and t+12, where 

number generator to generate a random number between 0 0<t<12. The marker sequence M at the end of the burst can 

and 9 and apply this number to select D. be used to resolve this modulo-12 ambiguity. 

In one embodiment, the maximum value of random delay In order to demodulate the payload part of the burst, the 

D is 9 symbols. Therefore, a user terminal can pick one of absolute timing information instead of the relative timing or 

the 10 random delays, {0,1,2, ... ,9} to send its CR burst 30 modulo-12 timing information is desired. As mentioned 

to the base station. In addition, the maximum allowable above, the marker sequence (M) at the end of the training 

roundtrip delay from the base station to the user terminal is sequence of FIG. 5 can be used for this task. The marker 

assumed to be 50 symbols. Therefore, a CR burst can be sequence M is set to be M=-x so mat one period of the core 

delayed as much as 59 symbols for a user terminal at the training sequence x augmented with M will be orthogonal to 

maximum possible range. This delay is measured with 35 two concatenated periods of x, i.e., [x, M] [x,x]"«0 in which 

respect to the CR burst of a user terminal that is at zero range x= t x i' ■ ■ • » x i2]> me 00115 sequence, forms a period of the 

and that did not delay its transmission. training sequence. 

The base station receives and resolves the CR burst as Since, in the present embodiment, the maximum allowed 

described below. The roundtrip delay is determined to be 4Q delay is 60 symbols, a roundtrip delay margin window with 

2(R/c)+D microseconds (since a symbol period is 2 micro- a duration of 5 periods of the core sequence is provided. This 

seconds where R km is the range of the user terminal. The leaves a duration of 4 periods of the core sequence for the 

base station, in the Configuration Message (CM) burst, then beamformer analysis window. Absolute timing is deter- 

instructs the user terrninal to advance its time by x a «(R/c+ mined by determining which of the four possible periods of 

D/2) microseconds in its response to the CR burst. The base 45 the core sequence is the last one received before the marker 

station is unaware of what, if any delay D the user terminal sequence. The beamformer analysis window may contain 

has applied to its CR burst. In response to this instruction, portions of five different core sequences, at least three of 

the user terminal, knowing its random delay D will subtract which are full periods of the core sequence. The first symbol 

out the random delay from the propagation delay and of the found sequence is either 12,24,36 or 48 symbols from 

advance its timing by x a -(D/2) microseconds. 5Q the first symbol of the marker sequence. Therefore the 

Hie random delay of D symbols is particularly valuable in absolute timing can be determined by testing only five 

a microcellular setting, or when user terminals cluster at hypotheses. In other words, given t, (0SK12), as the 

about the same range from the base station. If the range is modulo-12 delay estimate from the spatial processing 

about the same for two different user terminals, the CR search, the possible candidates for the absolute delay are: {%, 

bursts from the two user terminals will look identical at the 55 T+12 > t+24 > T+36 » 

base station and may not be resolved. With random delays, As can be seen from FIG. 5, for any delay less than 60 

the CR bursts will be spaced apart making them easier to symbols, the beamformer analysis window at the base 

distinguish. station will always contain at least 4 periods of the periodic 

sequence x. As a result, the first symbol, x lf of the first 

CR Burst Resolution 6Q repetition of the core sequence within the beamforming 

FIG. 5 shows a training sequence of a CR burst for one analysis window can be at most 12 symbols away from the 

embodiment of the invention. In FIG. 5, x510 is one of the border of the round trip delay margin window and the 

two core sequences {s 19 sj. x is the vector consisting of 12 beamformer analysis window. 

symbols x, to x 12 . In this application, boldface text is used This property of the repeating periodic training sequence 

to represent vectors and normal text is used to represent 65 can be taken advantage of when searching for the first 

scalars. The complete CR Burst training sequence is con- symbol, x Jf of the core sequence x. FIG. 6 shows a flow 

structed from the core periodic sequence x. Nine repetitions diagram for resolving the CR burst according to one 
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embodiment of the invention. The CR burst is received at the 
antenna array 602 and the measurements from the multiple 
antennas for the period that covers the beamformer analysis 
window are stored 604. Then the beamforming weights that 
will yield the minimum least squares (LS) error for each 
timing hypothesis in the search window of 12 symbols is 
determined. The least squares error results are used to find 
a coarse relative timing hypothesis 606. Finding this timing 
includes performing a coarse search for the symbol timing 
modulo- 12 and performing a follow-up search 608. 

A weight vector is then determined 610 using the relative 
timing. This weight vector is applied to the stored measure- 
ments to convert the measurements from each antenna 
channel into a single channel 612. Then, this single channel 
is analyzed to determine the fine timing 614 and the fre- 
quency offset 616 of the desired signal. The fine timing and 
frequency estimates are used to determine a new more 
accurate weight vector 618 based on the data in the beam- 
forming analysis window. This improves the quality of the 
weight vector and the new weight vector is applied to the 
stored measurements to convert the measurements from 
each antenna channel into a new single channel 620, further 
improving the quality of the received signal. The processing 
so far cannot perfectly determine the timing information, 
because it can only measure time delays modulo- 12 
symbols, the duration of the repeating core sequence. The 
timing estimate is in the range {0 . . . 12}, i.e., delays of x 
and x+12 symbols will ideally produce the same timing 
estimate at this point. Using the new weight vector, absolute 
timing is determined by finding the marker sequence 622. 
This resolves the modulo-12 timing ambiguity. 

When the absolute timing is known, the entire training 
sequence can be identified. Accordingly a new weight vector 
calculation can be made using the entire training sequence 
including the marker sequence 624. This redetermined 
weight vector is still more accurate because of the larger 
number of samples that can be used. It is applied to the 
stored measurements to again convert the measurements 
from each antenna channel into a single channel 626. With 
the new single channel, the CR burst is demodulated and 
read 628. 

This processing structure yields accurate beamformer 
weights that can be used to receive the desired signal without 
determining perfect timing for the signal. 

Coarse Search for Relative Timing 

After the CR burst has been received in each of the 
antennas of the array 602 and the signals have been stored 
604. The system can proceed to determine the relative 
timing based on the stored measurements 606. In one 
embodiment, this process is done using a covariancc matrix 
computation and a Cholesky decomposition. 

For measurements that lie in the beamforming analysis 
window assume that measurements are collected at 1.5 times 
the symbol rate. In this example, the beamforming analysis 
window is of duration AS symbols or 4 repetitions or periods 
of the core sequence. Next, assume that the variable core- 
Snap Point points to the first snapshot that is in the beam- 
former analysis window, and r(t) denotes the snapshot 
(measurement vector) collected at time t seconds after the 
start of the CR burst (relative to base station time). The 
covariance matrix R can be estimated as follows, where 
Thaw is lne time of the beamforming analysis window: 
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R = V riDr" (r), T B aw = U I coreSnapPolnr s t < coreSnapPoint+ 48} 

5 

After computing the covariance matrix R, a matrix L can 
be found such that 
R-LL* 

L is termed as the Cholesky factor of R. 

In one embodiment, the location of the first symbol of the 
core training sequence over the first 12 symbols of the 
beamforming search window is searched for using a least- 
squares processor with a search step of 2/3 symbols. This 
results in 18 possible locations of the first symbol, i.e. 18 
hypotheses and, accordingly, 18 least-squares (LS) error 
calculations. The vector coarseSearchGrid contains the 
delay hypothesis values for each 2/3 symbol increment (The 
unit used here is a symbol period): 
20 coarseSearchGrid-[0,2/3,4/3,2, . . . , 11 1/3] which con- 
tains the delay values for the 18 hypotheses. 

For each delay hypothesis k, (l^k§18), the following 
computations are performed: 

Calculate cross-correlation vector p k for each hypothesis: 

25 

30 In which r(t) is the received signal samples and d(t) is the 
samples of the desired signal. The desired signal samples are 
assumed to be available in a filter bank over-sampled at 24 
times the symbol rate. The value x k is the assumed delay for 
the k-th hypothesis, i.e., x^coarseSearchGrid (k). 

35 Back Substitution: The resulting vector p^is applied to 
solve for an interim vector x* for each hypothesis k, where 
L is the Cholesky factor as mentioned above: 

Lx*-p* 

Then the least-squares (LS) fit for each hypothesis is 
40 computed: 

/* = /(r*)= Ya M'+irf 

'* T BAW 

45 

Due to the periodicity of the desired signal within the 
beamforming analysis window, the cross correlation vector 
computation can be modified to: 

50 W/) + r(/ + 12) + r(f+24J + r(/ + 36)lrf»(f+r 4 ) 

in which T, is a period of the desired signal, i.e., 

55 ^-{tjcoreSnapPoint^tScoreSnapPomt+n} 

Therefore, the above calculation of the cross correlation 
vector p* costs one- fourth as many multiplies as it would if 
the core sequence were not a repeated periodic sequence 

60 within the training sequence. Finally, the hypothesis with the 
minimal LS error f* is picked. An illustration is presented 
with respect to FIG. 7. 

FIG. 7 shows the LS cost function on the vertical axis 702 
as a function of the time delay x for the 18 hypotheses on the 

65 horizontal axis 704. The horizontal axis is marked in units of 
2/3 of a symbol time in accord with the sampling rate of 1 .5 
times the symbol rate. The selected best coarse search time 
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delay estimate 706 is shown in FIG. 7 at t=4/3 symbols. 

Note that the delay is measured with respect to the nominal V iw i 2 *** 

time or modulo-12 time at which the CR burst is expected f{Jo) = h {du + T <> )] '^ x ° 
to arrive. This nominal time is named the coreSnapPoint. 
This simplified computation of the cross correlation vectors 

p* based on the periodicity of the training sequence within The beamforming operation results in a single channel 

the beamforming analysis window yields substantially the measurement that shall be named g(t) and defined as fol- 

same accuracy as a regular cross correlation computation. lows: 

Any mismatch of the training sequence in this simplified g(0" w «/' I (0 

approach makes only a minimal impact on the performance 10 The actual computations can be simplified because the 

of the beamformcrs due to the periodicity of the training time range of the beamforming operation does not need to be 

sequence. the whole burst. Note that the frequency estimation (search) 

has not yet been performed. This reduces the search to a 

rollow-up &earcn single dimension, namely time-delay. The frequency search 

Having made a coarse determination of the relative timing 15 can be avoided because of the duration of the beamforming 

606 a follow-up timing search 608 can be performed. This analysis window, the maximum frequency deviation 

search has one more LS error calculation. The result of this expected in the received signal and the target SNR level at 

calculation can be compared to the best LS error from the the output of the beamformer. When the beamformer analy- 

coarse search. sis window is short enough to make the frequency offset 

The LS error vector i k from the coarse search is treated as 20 negligible, the desired signal mismatch is Dot significant 

a circular buffer and concentration is on the midpoint of the enou 8 h t0 degrade performance. Furthermore, the inability 

range between the best location from the coarse search (that * m V™ c SNReven by using an MSE beamformer 

yielded the minimum LS error) and its neighbor that yielded ( that ^ows the perfect covariance matrix and steering 

the next best LS error " vector) to high levels (above 5 dB) makes the actual beam- 

. , ,' „„ ,. . , 25 former blind to imperfections in the desired signal (such as 

In the coarse search, the following 20 dimensional vectors fre<?llency ofbet) ^ presence of frequencv offset become s 

were iorme . serious at high target SNRs and ignoring large frequency 

updatedl^Values-[f 18 ,fj,f2,f3,f 4 , . . . , f 17 ,f J8 >fi] offsets causes the SNR performance to tend to saturate at 

and around 5 dB output SNR. However, this is still not an issue 

updatedCoarseSearchGrid-[-2/3,0,2/3,4/3,3 1/3, 4, . . . , 30 if the CR burst payload can withstand SNRs as low as 0 dB. 

11 1/3, 12] Furthermore, the single channel estimation of frequency and 

Assume that the hypothesis of f 3 (x-4/3) (of the current fine timing that follows will alleviate these problems using 

list) was selected by the coarse search procedure (as shown fewer computations, 

in FIG. 5) because f 3 <f2<f 4 . In the follow-up search . 

procedure, the LS error is computed for x-1, since this forms 35 hine J irmng bsliraalor 

the midpoint of the range between the delays corresponding The coarse search and the follow-up search provide a 

to f 2 (te2/3) and f 3 (t=4/3). This is shown graphically in FIG. timing delay with an effective rate of three times the symbol 

8. rate. The timing resolution can be improved to better 

FIG. 8 shows the LS cost function 802 as a function of the demodulate the payload. Severe intersymbol interference 

18 hypotheses 804 in which the time domain is again 40 can be expected at the worst case timing events (1/6 symbol 

sampled at a rate of 1.5 times the symbol rate. The selected period of timing error). To estimate fine timing 614, the 

best case coarse search time delay estimate 806 is x=4/3 Oerder-Meyr blind timing estimator can be used. This 

symbols. After the additional hypothesis test, the hypothesis estimator is blind in the sense that it does not assume any 

(among the 19 of them) which corresponds to the minimum knowledge of the transmitted symbols. Therefore it is insen- 

least -squares error is again picked. In this example, % 0 ~l is 45 sitive to frequency offset errors. The Oerder-Meyr estimator 

the delay value 808 that accomplishes this. can be used as follows: 

Beamforming Weights . ViTSt ! 6(0 to 3 times the symbol rate (from 1 5 

& 6 times the symbol rate) using an interpolation filter, and call 

After the follow-up search, determining the relative tim- 5Q the output g,<t). The duration of the measurements process 

ing delay x to within one third of a symbol period 608, the amounts to 48+12+82ol44 symbol periods, this is the guar- 

weight vector w 0 is recomputed to extract the desired signal anteed maximum number of symbols on the right hand side 

as follows: of the coarse liming estimate in the beamforming analysis 

Using the cross-correlation vector p 0 for the best bypoth- window as shown in FIG. 5. Alternatively, the symbols can 

esis: 55 be sampled at a three times rate originally. In such a case, 



only every other sample would be used in the coarse search 
Y (t)d (/ + r) and follow-up search. 

P ° ~ **tbaw r Next, the interpolated signal is passed through a memo- 

ryless nonlinearity to obtain 



And the corresponding interim vector x„: 



60 



Lx 0 «p 0 The absolute value nonlinearity generates tones at the 

The weights can be computed by solving the following multiples of the symbol rate. The nonlinearity is appro xi- 

back substitution: mated by fitting a polygon. For example, if z is a complex 

L // w c> -P 0 65 number with representation z-z r +jz„ then |z|-max(|z r j,|z l j)+ 

With the resulting LS error which is already determined to 0.34xmin flzj.lzj). The zero-crossing of the complex sine 

be f(x 0 ): wave at the symbol rate is determined by using a 3 point 
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DFT and this location is termed as a symbol transmission 
instant. Finally, the symbol transmission instant is found that 
is closest to the position identified by the follow-up search 
as the location for the first symbol of the periodic training 
sequence. 

The closed-form estimate obtained using the Oerder-Meyr 
timing estimator is termed as If this estimate is not in the 
range (0,12), it is mapped into this range by adding or 
subtracting 12 to it. 

Coarse Frequency Estimator 

After fine timing is acquired 614, the measurements are 
used in the analysis window together with the desired signal 
to compute the cross correlation function and determine the 
coarse frequency estimate 616 as the peak of the cross 
correlation function. For simplicity, this analysis is done at 
the 1 .5 times over sampling rate (instead of the 3 times rate 
used for the fine timing search), however many other sam- 
pling rates can be used. Specifically, for the m-th frequency 
hypothesis, the function freqMatch(m) is computed: 



freqMaichim) = £ *(0d*('+ r y ^)exp(-/2*/ wl /) 

The frequency candidate that maximizes freqMatch(m) is 
selected as the frequency estimate, fcoarse Frequency can- 
didates (f m *s) are in the set: 

{-2800,-2000,-1200,-400,400, 1200, 2000,2800} 
Alternatively, the frequency uncertainty can be regarded 
as much less and fewer hypothesis tests can be performed. 

More Accurate Beamforming 

Using the frequency and timing estimates t^, i coarse , an 
improved weight vector, w^,, can be determined 618 and 
applied to the found core sequences to obtain an estimate of 
the desired burst s(t): 

The more accurate weight vector is obtained by employ- 
ing d^+T^,) exp (J2jtf /iM4f t) as the desired signal in the LS 
error computations. However, if the corresponding LS error 
is larger than the follow up search based LS error f oJ the fine 
timing and frequency estimates are rejected and the esti- 
mates are set to % fifie ^ 0 and f cocrre .O, which results in 

To further improve the SNR results at high post-copy 
SNR levels, the search can be performed for a fine frequency 
estimate around the coarse frequency estimate at locations 
f coaw +[-800,0,800] using an LS fit to multichannel mea- 
surements for t eT^u,. Such a search can be included to 
improve the frequency estimation performance for the CR 
burst demodulator. 

Absolute Timing Recovery 

Prior to determining the absolute timing, as opposed to 
relative timing, the more accurate weight vector is applied in 
order to convert the measurements stored from each of the 
antennas into a single channel measurement 620. This 
conversion is still done using the data in the beamforming 
analysis window because only this data is certain to consist 
only of periods of the core sequence. An appropriate pro- 
cedure 622 for determining the absolute delay can be 
implemented as follows. Assume that the modulo- 12 time- 
delay estimate is t, and that g(t) denotes the output of a 
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beamformer (as described above g(t) is a single channel 
measurement) that extracts the desired CR burst. First, the 
candidate starting locations for the concatenated sequence 
[P9,M] (see FIG. 5) for each of the five delay hypotheses is 

5 determined. For example, if the absolute delay equals x, the 
sequence [P9,M] should start x+36 symbols away from the 
left border of the beamformer analysis window. For the 
hypothesis in which the absolute delay is assumed to be 
equal to r+48, the sequence [P9, M] should start x 84 

10 symbols away from the left border of the analysis window. 
More specifically, assume that T k denotes the begin time 
for the k-th hypothesis, i.e., T fc -x+Jf5+(12xk), for k-0, 
1,2,3,4. Let v k denote the samples of g(t) that correspond to 
the interval which starts at T k for a duration of 24 symbols 

15 (two periods). For each hypothesis, compute: 

M*)-M*-*]T- 

Pick the hypothesis that yields the maximum h(k). 

20 Assume that hypothesis k 0 is selected. Then, the absolute 
timing estimate is determined as T+(k„xl2). 

Due to the orthogonal nature of [x, m] and [x, x], there is 
a five-ary pulse -position modulation (PPM) demodulation 
problem in which the SNR changes depending on the true 

25 position of the desired vector [x, m]. For example, if the 
actual delay is less than a period, then the payload portion 
of the desired CR burst will appear as noise for the hypoth- 
esis tests 1,2,3,4. On the other hand, if the actual delay is, for 
example 20 symbols, then, the CR burst payload symbols 

30 will not contribute to the test of a hypothesis in which k=0 
as noise because of the orthogonality of [x, m] and [x, x], but 
they will contribute to the noise component in h(k) for 
k=2,3,4. With this delay condition, the beamforming analy- 
sis window receives the following 48 symbols: 

35 X 5> X 6» X 7» X 8» X 9> X 10,*11* X 12> X 1> X 2> X 3> X 4> • ■ • ^S>^6*^7»^8»^9» 
X l 0> X 1 1 12> X 1 > X 2> X 3 > X 4> 

In this stream, the first and last 12 symbols are written out. 
The spatial processing search will discover the first occur- 
rence of x,. at a delay of 8 symbols. The next step is to 
40 determine whether the absolute delay is 8,20,32,44,56 sym- 
bols. In order to determine this, the measurements following 
the beamforming analysis window are focused upon. With a 
delay of 20 symbols, the next 40 symbols will be: 

X 5> X 6j X 7> X 8> X 9» X 10» X J1> X 12' X 1> X 2» X 3> X 4> X 5> X 6> X 7> X 8> X 9> X 10» 
45 X J1 ,X 1 2,*"X 1 ,-Xj,— x 2 ,— X 3 ,— X 4 ,-X5,— X 6 , — x7,-x8,-Xp,- 

X 10»" X 11 F~Xl2idiA> d 3» d 4i d S 1 4«» d 7» d 8 

in which d, indicates the i-th payload symbol, which is 
unknown. As a result, v 0 =[x, x] and v^x, -x], which 
results in h(0)=0 and h(l)=24 2 . Accordingly, the 
50 desired signal does not contribute to the decision vari- 
able h(0) when the absolute delay is larger than 12. On 
the other hand, v^-x^-Xj, . . . ,-x 12> d^, . . . d 12 ] 
is not necessarily equal to zero. The first half of v 2 
equals the marker sequence, which yields the result: 

55 

A(2H-iMM* • • • 4 J*T- 
This differs significantly from: 

60 

h(2) is expected to be higher than h(3) on average because 
of the coherence term (-12) inside the magnitude operator 
and its presence creates a significant competition to the 
correct hypothesis decision value (h(0)) in the presence of 
65 measurement noise. Therefore, absolute delay errors that are 
equal to +12 should dominate other values of delay estima- 
tion errors. One way to eliminate the self-noise problem is 
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to put the marker sequence in the beginning of the burst and 
search for the concatenated sequence [-x, x]. In such a 
design, self noise will exist only due to the presence of the 
ramp-up symbols, but this will affect the noise component in 
only one of the h(k)'s. Other modifications to the ordering 
of the sequences are also possible to suit particular appli- 
cations. 

As discussed above, the remaining task is to find the 
marker sequence. This can be done based on a cross corre- 
lation operation with the 1.5 times over-sampled version of 
the sequence [x,-x] with v k corresponding to all of the five 
hypotheses locations indicating period slips of 0,1,2,3,4. 

Using the fine timing estimate x fing7 the time of the final 
beamformer output sample (that arrives at t c ) which is 
closest to coreSnapPoint+x^, is determined. Let x A denote 
the time difference between the chosen measurement at t c 
and the location of actual timing location: 

t^/c-coreSnapPoint-Tp,^ 

Remember that the beginning of the k-th hypothesis 
(k=0,l, 2,3,4) will be at times T^coreSnapPoint+x^. 
Therefore, the vector v k to be tested for the cross correlation 
for the k-th hypothesis will take the form: 

v*=[%+ 12x (3+Jfc)), £(f c +12x(3+fr)+2/3)L 
i(^+12x(3+Jt)+2x2/3), . . . 5(t c +12x(3+it>f 35x2/3) 

which is a vector with 36 entries. 

The v k 's are cross correlated with the desired signal vector 
d, which is formed from the knowledge of the training 
waveform and the timing difference x e : 

<^[j rf (ooreSnapPoint+T e + 1 2x3), ^coreSnapPo int+x e + 1 2x3+ 
2/3), . . . , jrf(corcSiiapPoint+T e +l 2x3+35x2/3)] 

in which s/t) is the training waveform for the desired signal. 

The magnitude squared cross correlations for each 
hypothesis is computed: 

A(*)-krfT 

Let k„ denote the index of the largest h(k). Then, the 
absolute timing is determined as: 

The values are in terras of symbol durations in the above 
equation with the period of the training sequence being 12 
symbols. 

Alternatively, symbol spaced sampling can be used in the 
absolute timing recovery. This can eliminate the mismatches 
created in the 1.5 times measurements that result in loss of 
orthogonality. With symbol spaced sampling, this operation 
will require 5 cross correlations of length 24 and comparing 
the squared-magnitudes of the cross correlations. 

After absolute timing recovery 622, the whole training 
sequence can be used to improve the timing and frequency 
estimates and these estimates can be used to perform one 
more stage of beamforming over the whole training 
sequence duration. This can increase the number of symbols 
used in spatial processing to 120 symbols from 48 and 
improve the performance at a cost of computational com- 
plexity The beamforming would involve first determining a 
still more accurate weight vector using the full training 
sequence 626 and then applying this weight vector to 
convert the measurements stored for each antenna channel 
into a single channel measurement 628. In doing so a new 
covariance matrix should be computed and therefore another 
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Cholesky factorization performed. The conversion can be 
applied to the entire CR burst so that the CR burst can be 
demodulated as a single channel 628. Alternatively, the 
absolute timing can be applied to determine the beginning 

5 and end of the CR burst and the existing weight vector and 
frequency estimates can be used to demodulate the signal. 
The choice will depend on the demands of the particular 
system and the condition of the channel. Having determined 
timing and a weight vector, this weight vector can be used 

10 to send traffic back to the remote from which the CR burst 
was received 630. 

General Matters 

The present invention can eliminate the need to acquire 
the absolute timing information in order to perform beam- 
forming. It exploits the periodicity of target training 
sequences in order to estimate the time delay encountered by 
the desired signal modulo-the period of the training 
sequence which is less than the total time delay ambiguity 

2Q range. The reception of uplink messages at the base station 
which are not perfectly synchronized is enhanced. 

Training based interference cancellation methods can be 
used to enhance the signal quality for such bursts. However, 
the training sequence should be aligned with the incoming 

25 signal, and therefore the resulting search process involves 
testing every possible delay with a predetermined resolution. 
As the delay ambiguity grows, this hypothesis testing pro- 
cess becomes increasingly more expensive with computa- 
tional resources. 

30 Using a periodic training sequence eases the efforts of 
hypothesis testing followed by a timing acquisition 
sequence to determine the perfect timing information. The 
periodic sequence decreases the timing ambiguity to within 
a period of the training sequence, which can be selected 

35 small so that resources are available to handle the hypothesis 
testing task. For example, if the expected delay is 50 
symbols, and if the period of the training sequence is 10 
symbols, with a resolution of a symbol, 10 hypotheses can 
be tested instead of 50, provided that the periodic training 

40 sequence can be observed in the search window, regardless 
of the delay. In addition, the search operation, which 
involves cross correlations can be simplified by the number 
of periods in the search window, if the search window 
contains 5 periods of the training sequence (period with 10), 

45 the cross correlations can be simplified by five-fold. 

Once the timing hypothesis is verified, only the number of 
periods of the repeating core sequence that has been skipped 
must be determined. To perform this more accurately, the 
multiple channel measurements (if available) can be con- 

50 verted to a single channel measurement using the search 
procedure results outlined above. After that a timing acqui- 
sition sequence that is at the beginning of the burst can be 
found, and to find the number of periods that were missed. 
So in the previous example, only periods {0,1,2,3,4} may 

55 have been missed in the beamforming analysis window. 
Accordingly, 5 searches can be performed and these are 
performed on single channel measurements rather than 
multichannel measurements. 
The invention eliminates the high cost of search on 

60 possibly multichannel data for all possible delays that may 
be encountered by the user terminal as a function of range. 
The burst structure and periodic training sequence makes it 
easier to search for the arrival time of the signal and the 
periodicity of the training further reduces the cost of com- 

65 putations. 

In the description above, for the purposes of explanation, 
numerous specific details are set forth in order to provide a 
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thorough understanding of the present invention. It will be determining timing for the burst based on the stored 

apparent, however, to one skilled in the art that the present measurements to which the determined weights have 

invention may be practiced without some of these specific been applied. 

details. In other instances, well-known structures and 2. The method of claim 1 wherein determining weights 

devices are shown in block diagram form. $ comprises selecting weights that will yield a minimum least 

The present invention includes various steps. The steps of squares error for at least a portion of the training sequence 

the present invention may be performed by hardware of the burst. 

components, such as those shown in FIGS. 1 and 2, or may 3. The method of claim 1 wherein the weights relate to the 

be embodied in machine-executable instructions, which may spatial directionality and interference of the received burst, 

be used to cause a general-purpose or special-purpose pro- Q 4. The method of claim 1 further comprising applying the 

cesser or logic circuits programmed with the instructions to determined timing to the burst and redetermining weights 

perform the steps. Alternatively, the steps may be performed that will yield a minimum least squares error for the burst 

by a combination of hardware and software. The steps have using the determined timing. 

been described as being performed by either the base station 5. The method of claim 4 further comprising applying the 

or the user terminal. However, any steps described as being redetermined weights to at least a portion of the burst and 

performed by the base station may be performed by the user 15 redetermining timing using the at least a portion of the burst 

terminal and vice versa. The invention is equally applicable to which the redetermined weights have been applied, 

to systems in which terminals communicate with each other c 6 ™ e m ^hod of claim 5 further comprising ^determining 

without either one being designated as a base station, a user offeet using the at least a portion of he burst to 

• i , , • t u<™k-, ct ,^„ which the redetermined weights have been applied, 

terminal, a remote terminal or a subscriber station. ^ ? ^ method Qf ^ 5 * ^ comprising applying the 

The present invention may be provided as a computer redeterm i ne d timing to the burst and further determining 

program product which may include a machine-readable we ights that wUl yield a minimum least squares error for the 

medium having stored thereon instructions which may be burst the redetermined timing, 

used to program a computer (or other electronic devices) to 8 ^ method of claim 5 wherein computing a covariance 

perform a process according to the present invention. The ^ ma trix comprises computing a covariance matrix only for a 

machine-readable medium may include, but is not limited to, portion of the known training sequence, 

floppy diskettes, optical disks, CD-ROMs, and magneto- 9. The method of claim 1 further comprising generating a 

optical disks, ROMs, RAMs, EPROMs, EEPROMs, magnet single representation of the received burst after applying the 

or optical cards, flash memory, or other type of media/ determined weights. 

machine-readable medium suitable for storing electronic 10. The method of claim 9 wherein determining timing 

instructions. Moreover, the present invention may also be further comprises calculating a cross correlation vector for 

downloaded as a computer program product, wherein the the plurality of different timing hypotheses for application to 

program may be transferred from a remote computer to a the least squares calculation. 

requesting computer by way of data signals embodied in a U. The method of claim 1 wherein determining the 

carrier wave or other propagation medium via a communi- weights comprises computing a covariance matrix against 

cation link (e.g., a modem or network connection). ^3 of claim t where in determining timing 

Importantly, while the present invention has been riscs determining a least squares error for a plurality of 

described in the context of a wireless internet data system for timing hypotheses for a portion of the known 

portable handsets, it can be applied to a wide variety of training sequence. 

different wireless systems in which data are exchanged. 4Q 13 ^ me thod of claim 1 wherein determining timing 

Such systems include voice, video, music, broadcast and comprises determining timing with respect to a repeated 

other types of data systems without external connections. corc sequence of the known training sequence and then 

The present invention can be applied to fixed remote termi- determining timing with respect to the entire received burst, 

nals as well as to low and high mobility terminals. Many of 14. The method of claim 13 further comprising redeter- 

the methods are described in their most basic form but steps 45 mining weights using the entire received burst after deter- 

can be added to or deleted from any of the methods and mining timing with respect to the entire received burst, 

information can be added or subtracted from any of the 15. The method of claim 1, wherein the burst is configured 

described messages without departing from the basic scope in accordance with at least one of a TDMA, a FDMA, a 

of the present invention. It will be apparent to those skilled CDMA and a FDD air interface standard, 

in the art that many further modifications and adaptations 5Q 16- A machine-readable medium having stored thereon 

can be made. The particular embodiments are not provided data representing sequences of instructions which, when 

to limit the invention but to illustrate it. The scope of the executed by a machine, cause the machine to perform 

present invention is not to be determined by the specific operations comprising: 

examples provided above but only by the claims below. receiving a burst from a terminal at a set of spatial 

What is claimed is: 5S diversity antennas, the burst including a known training 

1, A method comprising: sequence; 

receiving a burst from a terminal at a set of spatial storing measurements of the received burst from each 

diversity antennas, the burst including a known training spatial diversity antenna corresponding to receive times 

sequence; within a selected analysis window; 

storing measurements of the received burst from each 60 determining weights for each diversity antenna for at least 

spatial diversity antenna corresponding to receive times a portion of the training sequence of the burst before 

within a selected analysis window; timing for the burst has been determined; 

determining weights for each diversity antenna for at least applying the determined weights to the stored measure- 

a portion of the training sequence of the burst before ments; and 

timing for the burst has been determined; 65 determining timing for the burst based on the stored 

applying the determined weights to the stored measure- measurements to which the determined weights have 

ments; and been applied. 
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17. The medium of claim 16 wherein the instructions 
causing the machine to perform operations comprising 
determining weights further comprise instructions for select- 
ing weights that will yield a minimum least squares error for 

at least a portion of the training sequence of the burst. S 

18. The medium of claim 16 further comprising instruc- 
tions causing the machine to perform operations comprising 
applying the determined timing to the burst and redetermin- 
ing weights that will yield a minimum least squares error for 
the burst using the determined timing. 10 

19. The medium of claim 18 further comprising instruc- 
tions causing the machine to perform operations comprising 
applying the redetermined weights to at least a portion of the 
burst and redetermining timing using the at least a portion of 
the burst to which the redetermined weights have been 15 
applied. 

20. The medium of claim 18 further comprising instruc- 
tions causing the machine to perform operations comprising 
determining frequency offset using the at least a portion of 
the burst to which the redetermined weights have been 20 
applied. 

21. The medium of claim 18 further comprising instruc- 
tions causing the machine to perform operations comprising 
applying the redetermined timing to the burst and further 
determining weights that will yield a minimum least squares 25 
error for the burst using the redetermined timing. 

22. The medium of claim 16 wherein the instructions 
causing the machine to perform operations comprising 
determining the weights further comprise instructions to 
compute a covariance matrix against the expected signal for 30 
a portion of the known training sequence. 

23. The medium of claim 16 wherein the instructions 
causing the machine to perform operations comprising 
determining timing further comprise instructions for deter- 
mining a least squares error for a plurality of different timing 35 
hypotheses for a portion of the known training sequence. 

24. The medium of claim 16 wherein the instructions 
causing the machine to perform operations comprising 
determining timing further comprise instructions for deter- 
mining timing with respect to a repeated core sequence of 
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the known training sequence and then determining timing 
with respect to the entire received burst. 

25. The medium of claim 16, wherein the burst is con- 
figured in accordance with at least one of a TDMA, a 
FDMA, a CDMA, and a FDD air interface standard. 

26. An apparatus comprising: 

a set of spatial diversity antennas at a terminal to receive 
a burst, the burst including a known training sequence; 
a memory to store measurements of the received burst 
from each spatial diversity antenna corresponding to 
receive times within a selected analysis window; 
a processor to determine weights for each diversity 
antenna for at least a portion of the training sequence of 
the burst before timing for the burst has been deter- 
mined to apply the determined weights to the stored 
measurements, and to determine timing for the burst 
based on the stored measurements to which the deter- 
mined weights have been applied. 

27. The apparatus of claim 26 wherein the processor 
further applies the determined timing to the burst and 
redetermines weights that will yield a minimum least 
squares error for the burst using the determined timing. 

28. The apparatus of claim 27 wherein the processor 
further applies the redetermined weights to at least a portion 
of the burst and redetermines timing using the at least a 
portion of the burst to whichfihe redetermined weights have 
been applied. 

29. The apparatus of claim 27 wherein the processor 
further determines frequency offset using the at least a 
portion of the burst to which the redetermined weights have 
been applied. 

30. The apparatus of claim 27 twherein the processor 
determines timing with respect to a repeated core sequence 
of the known training sequence and then determines timing 
with respect to the entire received burst. 

31. The apparatus of claim 26, wherein the terminal is 
within at least one of a TDMA, a FDMA, a CDMA and a 
FDD radio communications system. 

***** 
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